Topological Phases in Graphitic Cones 
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The electronic structure of graphitic cones exhibits distinctive topological features associated 
with the apical disclinations. Aharonov-Bohm magnetoconductance oscillations (period $0) are 
completely absent in rings fabricated from cones with a single pentagonal disclination. Close to 
the apex, the local density of states changes qualitatively, either developing a cusp which drops to 
zero at the Fermi energy, or forming a region of nonzero density across Ef, a local metalization of 
graphene. 
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If one or more sectors are excised from a single layer of 
graphite (hereafter, graphene) as illustrated in Figure 
and the remainder is joined seamlessly, a cone results. 
Researchers discovered how to produce such graphitic 
cones in 19970]. We show that the topological nature of 
the apical defect profoundly modifies the low-energy elec- 
tronic structure, with remarkable consequences for non- 
local transport properties in this new class of nanoscale 
carbon. Rings made from cones exhibit an anomalous 
Aharonov-Bohm effect. In some cones, the fundamental 
period is completely absent. Also, the Fermi-level density 
of states near the apex is enhanced and shows distinctive 
energy dependence. 

We classify a cone by its opening angle, or equivalently 
the number of sectors removed. All possible angles 
have been observed experimentally; we concentrate on 
= 1 and = 2^ which reveal the essential physics. 
Since pentagonal defects are expected to predominate in 
real cones, is also typically equal to the number of 
defects, which must be tightly clustered to produce a 
conical shape. 

The gap of a semiconducting tube arises from 
frustration!!, ^, ^ in the phase of the graphitic Fermi- 
level states as they pass around the tube circumference. 
In this sense, the semiconducting nanotube is a transla- 
tional phase defect. A cone forms a rotational defect 
which also frustrates electronic phase. In nanotubes, 
the frustration affects states at the two graphene Fermi 
points in precisely the same way. In contrast, the topo- 
logical phase in a cone of odd entwines the two Fermi 
points. Although most of a large cone is just a gently 
curved graphene surface, the global electronic properties 
are significantly disturbed. 

Before presenting the new results, we briefly recall the 
continuum theory]^, ||, |^ of graphene, as it forms the 
framework for this understanding. We require a flexible 
real-space continuum formulation to handle structures 
which are locally graphitic but contain non-contractible 
closed loops, i.e., those surrounding the waist of a nan- 
otube or the apex of a cone. We think of the entire cone 
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FIG. 1: Constructing a cone from a graphene sheet. Cut 
out sector /3ay and join a/? to 07 to produce a {n^ = 1) cone 
with an apical pentagon. Excise a second sector and join afl 
to aS to produce a {n^ = 2) cone with a smaller opening 
angle. The resulting four-membered ring is very unstable; 
actual cones of that opening angle almost surely contain two 
nearby pentagonal defects instead. 



as having a perfect graphene-like structure locally, by 
considering a pentagonal defect not as a five-membered 
ring, but as a hole in a graphene structure. In the contin- 
uum picture this makes sense and the hole can be shrunk 
to a point. We can then represent defects as pure gauge 
fields. Our continuum description is ideally suited to 
bringing out topological aspects, so long as one exercises 
care in handling phenomena on atomic length scales. 

Graphene's Fermi surface comprises only two points. 
Figure || depicts the Fermi point states in a tight-binding 
language. The states are labelled by the direction of a 
wave-vector K and a sublattice index A or B. The figure 
illustrates \KA) and \KB), where K points up the page. 
The Fermi level space is four-dimensional, and we choose 
\K±A) and \K±B) as our basis, with K- = Un- 
der rotation counterclockwise about any hexagon center 
by 120°, implemented by i?6^, these states pick up simple 
phases: Re^\KA) = r]\KA), and Re^\KB) = r]\KB), 
where 77 — exp(27rz/3). 
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FIG. 2: Tight-binding wavefunctions of graphene Fermi- 
level states. Pick a direction K pointing from the center of 
a hexagon to the middle of one side, and assign amplitudes 
as shown, where 77 = exp(27ri/3). If the nonzero amplitudes 
are on the right-hand atoms of bonds cut by K, as in (a), 
the state is \KA), if on the left-hand atoms, as in (b), \KB). 
There are only two distinct K. 



We avoid momentum space notions because they are 
not available in systems such as cones, which lack trans- 
lational invariance. K is only a label for the orienta- 
tion of the pattern of amplitudes. However, whether in 
a flat graphene sheet, a nanotube, a cone, or any other 
graphitic structure, the states near the Fermi level must 
look like the four we have constructed, at least locally. 
Globally, low-energy states consist of slowly-varying en- 
velope functions multiplying these four patterns. We col- 
lect the envelopes into a four-component wavefunction 
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where the ± subscripts refer to and K . This is 
very similar to an effective mass description. We discuss 
obstructions to a globally consistent choice of later. 
If that choice is successfully made, A and B present 
no further difficulty. 

To complete the picture, we need a Hamiltonian. To 
avoid complicating distractions, we first consider simple 
graphene wherein crystal momentum is a good quantum 
number and we can work with just the components. 
In that case, the Hamiltonian is[^ ^, ^ 
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The cr's are Pauli spin matrices operating in the A/B 
indices; A is 'up,' B is 'down,' and the x axis is along 
K^. For the components of the wavefunction, the 
appropriate frame [x and y axes) is rotated through 180° 
relative to that for However, by an appropriate 

change of phases (±7r/2), the Kj^ frame may be used for 
all components. The frame need not be locked to the un- 
derlying lattice, however, a fact which is very useful on 
nontrivial surfaces such as cones. If the coordinate frame 
{e^,, e^} in Eq. (||) is rotated by 9 counterclockwise, and 
the wavefunction simultaneously transformed by multi- 
plication by exp (i0(T3/2), the form of the Hamiltonian is 
preserved. This invariance is made local, by introducing 
a gauge field, replacing d by d — ia, where a = —0-386/2. 



We imagine taking the sheet in Fig. ^with wedge Paj 
missing and wrapping it up to form a cone. The edge 
is then labelled by azimuthal angle (j) — and aP by 
(j) = 27r. What are the appropriate boundary conditions? 
Using the four Fermi point patterns described above, one 
can see that, up to phases, \K+A;(p ~ 0} is the same 
thing as \K-A; (j> = 2tt), and \K-B; (p — 0) the same as 
\K+B; 4> ~ Stt). There is no globally coherent distinction 
between and K^, a remarkable conclusion. In a 
sense, just as a Mobius strip has only one side, cones 
with odd have only one Fermi point! 

Before dealing with the awkward boundary conditions, 
notice that the local frames {ex^ey} are discontinuous 
across the seam. Using the local rotation invariance men- 
tioned earlier, we orient our frames so that — 
points tangentially, and Sy = — toward the apex. This 
change introduces a gauge field (now for general n^^) 
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where £ = [1 — n^/&\2TTr is the circumference of the cone 
at distance r from the apex. The adjustment of the 
frames makes an additional contribution to the bound- 
ary conditions. The various ingredients combine to give 
^'(^ = 27r) = T{n^/)'^{(j) = 0), where the operator T (a 
holonomy) is 



r(nj = (-l)exp 
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T2 is the standard (2 ^ y) Pauli matrix acting on the 
Fermi point indices ±. For odd , the exponential factor 
is eliminated by a singular gauge transformation, d — 
ia ^ d — ia — if3 with 

Recall that in solving the ordinary Schrodinger equation 
on a flux-threaded ring, the vector potential can be elim- 
inated by imposing a discontinuous boundary condition 
i/j{9 = 2tt) = exp[ie/ A ■ dl]tp{0). Here we perform es- 
sentially the reverse procedure, but our gauge field is 
proportional to a Pauli matrix. For even , the two 
Fermi points are not mixed and the exponential factor is 
just —1. So the gauge transformation need not involve 
T-matrices, but merely introduces half a flux quantum of 
fake magnetic flux. We deal with the first factor of — 1 
in Eq. (||) by using anti-periodic boundary conditions. 

With the rotational symmetry more manifest, we can 
now make a partial-wave decomposition of a general 
spinor as 
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and a similar decomposition of the Hamiltonian as H = 
^ h^^K The total angular momentum takes on all half- 



integer values, j 



^3/2,-1/2,1/2, 



We work 
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with a finite cone, so that the radial wavefunction x*^-') 
is in the Hilbert space of functions on (0, R] which are 
square integrable with respect to r dr, and the radial 
Hamiltonian h^^'' is of the form 



hi, = v„ 
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dr 
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The value of v depends on j and the type of cone involved: 
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FIG. 3: Schematic densities of states for a small patch near 
the apex of a cone, according to Eq. (^). 




n„ = 1; 
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We include the flat sheet (n^^ — 0, no disclination) as a 
'control,' and also introduce a magnetic flux <I> through 
the disclination in order to study magnetoconductance 
(<i>o = h/e is the normal flux quantum). Notice T2 in the 
n^j = I expression — energy eigenstates are superposi- 
tions of the two K. 

Comparing the eigenvalue equation for h„ to standard 
recursion relations for cylinder functions |l^ reveals 
the solutions to be ordinary Bessel functions: 
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The indicated restrictions on v involve a subtlety, ft is 
insufficient simply to ask that x a-nd h^x be square inte- 
grable. That requirement says that Xv.i is acceptable for 
— 1/2 < and Xu,2 for v < 1/2. Also requiring that the 
radial Hamiltonian h^, be self-adjoint, as it must, forbids 
both Xu,i and Xiy.2 in the domain of h^, at the same time, 
except for the special case v — 0. 

Now we turn to observable consequences. The energy 
density of states near the apex of a cone shows a remark- 
able dependence on the opening angle of the cone. In 
fact, the strictly local density of states diverges as r ^ 
for most cones, so we investigate instead the total den- 
sity of states on a patch < r < S for small 5. This 
quantity is more relevant anyway for comparison to both 
experiments and tight-binding computations. The result 
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as illustrated schematically in figure ^ 

Eq. (^) is derived as follows. From the large- a; asymp- 



totic formula ~ {2/'kx) 



1/2 



cos X ■ 



1/2 



iTrk/2Ry'^ for our 



deduce a normalization factor of c 
nificantly to the density of states near r = are j = ±1/2 



functions. The only values of j which contribute sig- 
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0,1 and j = 1/2 for ' 



2, and the leading or- 



der behavior of J„ is used to integrate over a small disk. 
Finally, the spacing of k values is Ak — tt/R {2tt/R for 
j = 1/2 and = 2). Thus, the total density of states 
in the ^-disk contributed by states associated with Bessel 
function order n is 
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Inserting appropriate values for n, we get Eq. (|9|). 

Strikingly, the low-energy density of states for the n^^ — 
1 cone has a cusp at Ep, and the cone with two sectors 
missing has a non-zero density of states near the Fermi 
level. An actual cone with the = 2 opening angle 
probably actually has a pair of = 1 defects. However, 
if the two pentagons are very near each other, we expect 
this conclusion to remain valid. These remarkable results 
should be observable in STM spectroscopy, so long as 
they are distinguished from the local effect of bond strain 
in the pentagonal defect (s). 

Although the continuum description lacks atomistic 
detail, the qualitative conclusions agree well with an ear- 
lier atomistic computationpl of (apical) local densities 
of states using a one-orbital tight-binding model for sim- 
ilar structures. Direct comparison with our Eq. can 
be made, bearing in mind the low-energy restriction on 
our calculations. In references ^ it was suggested that 
the states contributing to the non-zero density of states 
at Ej^ for rifj = 2 are (power-law) localized. Our cal- 
culations show that they are extended states which are 
enhanced in the vicinity of the apex. 

A mesoscopic normal metal ring in a perpendicular 
magnetic field exhibits oscillations as the flux through the 
ring is varied JT^. The longest period of oscillation is one 
normal flux quantum <I>o — h/e, corresponding essentially 
to the Aharonov-Bohm effect, and first observed in this 
context in the mid-1980's||l^. Manufacturing such a ring 
from a graphitic cone seems possible. The apex could be 
etched away with acid, or cut off with an STM. Such 
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manipulation does not change because it measures a 
topological property — the number of missing sectors, or 
equivalently, the opening angle. 

The magnetoconductance of such conical rings also 
shows profound sensitivity to the opening angle. Al- 
though we could argue directly from the holonomies, it 
is perhaps simpler to appeal to Eq. (^). These expres- 
sions are still valid for a ring, even though the subse- 
quent analysis was specific to a cone with only the apical 
point removed. For n^^ — and n^^ — 2, the Hamilto- 
nian is SU{2)k invariant (the SU{2) group acts on the 
K± indices), so that the two Fermi points produce two 
independent branches of excitations which respond iden- 
tically to a magnetic flux. Each branch exhibits its own 
magnetoconductance, which has period $o, just as in an 
ordinary metal, and the phases are identical. The total 
response therefore also has fundamental period $o- 

In ann^ = 1 ring, on the other hand, the $o component 
of the oscillations is extinguished. From Eq. (^) , we see 
that the T2 = +1 (r2 = — 1) branch behaves as though it 
were subjected to a flux of $ = <I>o/4 = -$o/4). This 
relative shift of the magnetoconductance curves causes 
cancellation of the $0 periodic component. Cones with 
intact apices will also show novel magnetic phenomena, 
including the prospect for field-tuneable radial charge 
density waves. 

Our neglect of disorder and inelastic scattering in this 
discussion should be permissible up to the micron length 
scale. The elastic scattering length, ip, in single- wall 
nanotubespUl is believed to range up to many microns 
and there is direct evidence that the phase coherence 
length {£^) at room temperature is also that longp^. A 
graphitic ring would be expected to have a somewhat re- 
duced £p, due to rough edges, but is likely comparable 
to that in the tube. 

In conclusion, we demonstrate that graphene cones, 
which have been experimentally produced, but not yet 
adequately studied, comprise a new class of nanoscale 
carbon wherein phase frustration induces profound mod- 
ulations of the low-energy electronic properties. These 
rotational phase defects (induced by topological lattice 
defects) leave distinctive and non-trivial local and global 
imprints on the electronic structure of graphitic cones. 
We predict two specific phenomena flowing from this ob- 
servation: apical enhancement of density of states and 



an anomalous magnetoconductance. 
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